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Abstract 

discuss  the  lattice  dynamics  of  argon,  krypton  and  xenon  overlayers 
on  the  Ag(111)  surface.  We  consider  monolayer,  bilayer,  trilayer,  and 
twenty-five  layer  films  of  each  of  these  adsorbates.  Data  are  also  presented 
on  the  dispersion  relation  of  selected  branches  of  the  phonon  spectra  of 
these  overlayers.  The  data  have  been  obtained  by  the  method  of  angle-resolved 
Inelastic  He  scattering.  Several  models  of  the  lattice  dynamics  are  compared 
with  the  data.  It  Is  concluded  that  the  gas  phase  potentials  proposed  by 
Barker  and  co-workers  prove  suitable  for  a  description  of  lateral  Inter¬ 
actions  between  the  adsorbates,  within  the  accuracy  of  the  available  data, 
provided  that  the  phonon  spectra  are  calculated  for  a  lattice  with  the 
experimentally  determined  lattice  constants.  '\  o  \  P  o  T  '  J  '* 
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I.  Introduction 

Oaring  the  past  two  or  throe  years ,  the  dispersion  relations 
of  surface  phonons  and  surfaee  resonance  nodes  have  been  studied 
experimentally  on  several  crystals.  Three  methods  have  been 
used.  The  first  is  the  inelastic  scattering  of  highly  aonoenes- 
getic  beams  of  helium  atoms  from  surfaces.1.  The  second ,  neutron 
scattering/  has  been  done  on  adsorbed  monolayers/  but  this  method 
only  works  for  high  surface  area  substrates  (for  example/ 
exfoliated  graphite) ,  and  only  for  monolayers  composed  of 
mstsrials  with  good  coherent  neutron  scattering  cross  sections. 

Of  the  rare  gases,  only  ir3<  falls  into  the  latter  category2. 

The  third  is  the  inelastic  scattering  of  low  energy  electron 
beams.3  Oil  three  techniques  are  very  similar  from  the 
conceptual  point  of  view,  and  provide  as  with  realizations  of  a 
surface  analogue  of  the  inelastic  neutron  scattering  experiments  . 
that  have  played  such  a  crucial  role  in  elucidating  the  nature  of 
elementary  excitations  in  bulk  crystals,  we  now  have  detailed 
data  in  hand  on  a  number  of  alkali  halide  surfaces/1  those  of 

i 

noble  metals z1  and  both  clean3  tad  adsorbate  covered*' 5  surfaces 
of  transition  metals. 

This  paper  discusses  data  on  the  surface  phonon  dispersion 
curves  associated  with  systems  distinctly  different  in  nature 
from  those  just  mentioned.  These  are  ordered  over layers  of  the 
heavy  tare  gases  (argon,  krypton,  xenon)  on  the  Ag(lll)  surface. 
The  emphasis  here  is  on  the  vibrational  motions  of  the  rare  gas 
adsorbates.  The  rare  gases  are  physisorbed,  rather  then  chemi¬ 
sorbed  on  this  surface.  The  physisorption  potential  well  is  both 


i 
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shallow  sad  broad#  with  the  consequence  that  tha  restoring  forca 
for  vertical  motion  of  tha  adlayer  is  vary  weak#  and  tha 
aaaoeiatad  vibrational  frequency  ia  low.  At  tha  sona  boundary  of 
tha  two  dimensional  Brillouin  sona#  tha  phonon  fraquaneiaa  of 
thaaa  adlayars  lia  wall  balow  thoaa  of  tha  substrata#  in 
constrast  to  thoaa  aaaoeiatad  with  light  chemisorbed  atoms#  which 
tand  to  lia  wall  above  tha  substrata  phonon  bands  throughout  tha 
two  diaanaional  Brillouin  zona.*  Haliua  baaas  of far  a  unique 
a sans  of  probing  such  vary  soft  surface  phonons#  whose 
fraquaneiaa  lia  balow  30  ea~^  throughout  tha  Brillouin  sona  for 
tha  systaaa  considered  hare.  Tha  energy  resolution  available 
from  tha  inalastie  electron  scattering  method  is  as  yet 
insuffieiant  to  enable  resolution  of  these  nodes. 

Extensive  studies  of  such  rare  gas  overlayers  on  Ag(lll)  have 
bean  reported  by  Brueh#  Webb#  and  co-workers . 7 * 8  Tha  physical 
picture  which  aaargas  is  tha  following.  From  tha  perspective  of 
tha  rare  gas  atoa#  tha  Ag(lll)  surface  looks  perfectly  smooth; 
tha  substrata  provides  only  tha  physisorption  potential,  7 (z), 
which  binds  tha  rare  gas  atoms  to  tha  surface#  and  there  is  no 

evidence  for  any  influence  of  tha  apparently  vary  weak  modulation 

•  • 

of  this  potential  in  tha  two  directions  parallel  to  tha  surface. 
Thus#  tha  adsorbates  form  hexagonal  overlayers  incommensurate 
with  tha  substrata#  apparently  with  lattice  constant  controlled 
only  principally  by  lateral  interactions  between  the  rare  gas 
atoms.  The  lattice  constants  measured  are  quite  close  .  though 
not  precisely  equal,  to  those  calculated  through  the.  me  of  rare  gas 
interaction  potentials  deduced  from  gas  phase  data. 
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This  paper  reports  data  on  the  phonon  dispersion  curves  of 
argon,  krypton  and  xenon  monolayers,  bi layers  and  tri layers  on 
the  Ag(lll)  surface.  In  addition,  we  report  data  on  films  formed 
from  approximately  twenty  five  layers  of  adsorbate;  to  the  helium 
atom,  these  apear  very  similar  to  semi-infinite  crystals.  We 
also  present  calculations  of  the  dispersion  curves  of  such 
structures  for  a  range  of  models  of  various  levels  of  sophisti¬ 
cation.  Throughout  the  paper  we  assume,  following  Webb  and  co¬ 
workers,  that  the  substrate  is  passive,  and  its  only  role  is  to 
provide  the  physisorption  potential  V(z).  We  then  explore  the 
sensitivity  of  the  calculated  phonon  dispersion  curves  to  varia¬ 
tions  in  the  lattice  dynamical  models,  which  range  from  the 
nearest  neighbor  central  force  models  applied  to  unrelaxed 
multilayer  structures,  to  calculations  on  fully  relaxed  multi¬ 
layer  structures,  with  lateral  interactions  modeled  through  use 
of  current  rare  gas  pair  potentials  extracted  from  gas  phase 
data.  The  results  are  compared  with  the  data,  to  test  the 
adequacy  of  such  pair  potentials  for  describing  lateral  inter¬ 
actions  in  the  adsorbed  state. 

We  should  comment  on  the  philosophy  of  the  theoretical 
analysis  contained  in  the  present  paper.  For  the  rare  gas 
adsorbates  considered  here,  excellent  gas  phase  pair  potentials 
are  available,  as  indicated  in  the  preceding  paragraph.  Of 
course,  these  pair  potentials  are  necessarily  modified  somewhat 
in  the  course  of  adsorption.  It  is  very  difficult  to  estimate 
these  changes  in  a  quantitatively  reliable  manner.  Also,  new 
sources  of  lateral  interaction  can  occur,  such  as  those 
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associated  with  adsorption  induced  dynamic  electric  dipole 
moments  (whose  strength  is  unrelated  to  the  static  dipole  moments 
deduced  from  work  function  changes).  Xn  addition,  when 
multilayer  films  are  considered,  the  physisorption  potential 
experienced  by  adsorbates  in  the  second  and  third  layer  differs 
from  that  of  adsorbates  in  the  innermost  layer,  because  of 
screening  of  the  van  der  Waals  interaction  by  the  inner  layers. 

This  is  also  difficult  to  quantify  reliably,  at  the  microscopic 
distances  of  interest  here.  Finally,  the  influence  of  motions  of 
substrate  atoms  can  induce  frequency  shifts  of  the  adsorbate 
phonons.  A  number  of  authors  have  explored  coupling  of  adsorbate 
and  substrate  atom  motions,9  but  the  calculations  are  not  applicable 
to  the  systems  of  present  interest  unfortunately.  Elsewhere,  within 
a  simple  model,  we  have  discussed  this  question  explicitly10  for  the 
present  systems.  There  are  indeed  frequency  shifts  produced 
comparable  to  those  associated  with  the  above  mentioned  corrections. 

We  thus  proceed  by  inquiring  whether  these  modifications  are 
crucial  by  simply  using  the  unmodified  gas  phase  pair  potentials  to 
model  the  lateral  interactions,  and  combining  this  with  the  Cole- 
Vidali  form  of  the  physisorption  potential  appropriate  to  an 
isolated  adatom.  We  find,  in  fact,  that  the  picture  just  described 
works  very  well  indeed,  when  the  experimentally  measured  lattice 
constant  is  employed  in  the  calculations.  We  also  explore  the 
phonon  spectra  generated  from  some  simpler  models,  such  as  a  nearest 
neighbor  central  force  model,  and  also  use  of  the  Lennard-Jones 
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potential  to  model  lateral  interactions,  and  we  find  these  clearly 
inadequate. 

Section  II  of  this  paper  discusses  the  experimental  procedure 
and  the  data,  and  Section  III  discusses  our  theoretical  calcula¬ 
tions,  and  comparisons  between  theory  and  experiment.  Section  IV  is 
devoted  to  conclusions  which  may  be  drawn  from  this  work. 


1 1 .  Experimental  Procedure 


The  surface  phonon  dispersion  curves  discussed  in  this  paper 
wer e  obtainod  using  angle  resolved  inelastic  Ho  scat taring  aa  tho 
dynamical  psobo.  Tho  actual  time-of- flight  (TOF)  measurements 
were  earsiod  out  ia  aa  ultra-high  vacuum  acattoriag  apparatus 
having  a  base  proaauro  o£  -  2  *  10*10  tors.  Ia  this  apparatus 
tho  incident  and  final  acattoriag  aagloa  can  bo  independently 
varied/  thus  allowing  tho  Inelastic  acattoriag  TO F  spectra  to  bo 
collected  across  tho  satire  Bsillouia  rone  using  a  fired  angle  of 
incidence.  This  is  an  important  feature  of  the  present  experi¬ 
mental  arrangement  for  accurately  determining  inelastic  scatter¬ 
ing  probabilities.  Detection  of  the  scattered  He  atoms  was 
accomplished  using  a  differentially  pumped  gnadrupole  ossa 
spectrometer#  for  which  the  crystal  to  ionise r  distance  *«•  14.45 
cm.  The  angular  resolution  of  the  detector  was  —  0.67°.  A  more 
detailed  description  of  the  experimental  apparatus  and  procedures 
will  be  presented  elsewhere  [n]. 

. The  substrate  for  these  experiments  was  a  single  crystal  of 
Ag#  cut  and  polished  to  within  0.5°  of  the  (111)  face,  as  checked 
by  Laue  X-ray  back-reflection.  The  crystal  surface  was  cleaned 
using  Ar*  ion  bombardment.  Cleanliness  was  verified  using  Auger 
spectroscopy .  After  annealing,  the  surface  order  and  orientation 
were  checked  using  H2  diffraction.  Surface  coherence  was  also 
checked  by  analyzing  the  specular  reflection  angular  profile 
using  a  63  meV  He  beam.  The  coherence  length  was  found  to  be 
on  the  order  of  100  A  after  annealing  at  750°  K.  Temperature  was 
monitored  bv  a  Chromel-Alumel  thermocouple  spot-welded  directly 


7 


to  the  back  of  the  crystal.  The  crystal  temperature  was  bald 
at  2S  t  1°X  for  Xr  and  Xa  measurement s,  and  at  21  ±  1°K  for  Ac 
measurements . 

A  mw  liquid  nitroqan  coo lad  baas  sourca  was  usad  to  produca 
a  vary  stabla  13  amV  Ha  bams  [n].  Tha  velocity  dispersion,  av/t, 
at  this  energy  was  -*  It.  Tha  incidant  basis  was  cBachanically 
choppad  by  a  rotating  slottad  disk  having  a  shuttar  function  of 
3  iisac  FWHM.  Tha  baas  was  collimated  to  an  angular  divarganca  of 
0*3°. 

Tha  thin  film  growth  procedures  davalopad  during  tha  aarly 
stag as  of  thasa  experiments  wars  guidad  by  tha  larga  quantity  of 
structural  and  thermodynamic  data  that  had  baan  praviously 
coapilad  for  cars  gas  ovarlayars  physisorbad  on  Ag(lll)  [7,3]. 

In  particular,  rara  gasas  fora  asiauthally  aligned  but  trans- 
lationally  incoaaansurata  structures  on  this  surface.  Ovarlayars 
ware  grown  using  a  differentially  pumped  dosing  beam  inclined  at 
154  to  tha  scattering  plana.  X  full  recipe  for  tha  rather 
complicated  process  of  growing  good  ovarlayers  will  be  presented 
in  a  subsequent  paper  [11].  However,  it  is  important  to  note  that 
our  ovarlayer  growth  procedure  is  not  based  upon  thin  fU.m  growth 
occurring  at  equilibrium.  In  brief,  the  Ar,  Kr,  and  Xe  mono- 
layers  ware  grown  by  controlling  the  dosing  beam  pressure  and  the 
substrate  temperature.  Multilayer  films  could  then  be  grown  on 
tha  initially  deposited  and  carefully  anealed  monolayers  using 
computer  controlled  dosing  times  and  a  multistep  annealing 
procedure.  The  long-range  order  and  lateral  spacings  ware 
checked  by  elastic  Ha  diffraction;  examples  are  shown  in  Figure  l 
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for  1,  2,  3,  and  25  layer  (111)  oriented  films  of  Ar.  Analysis  of  the 
specular  peak  widths  shown  in  Figure  1  indicates  that  all  of  the  rare  gas 
films  have  coherence  lengths  that  are  only  slightly  smaller  than  the  Ag 
substrate,  *  100  X.  The  very  low  diffuse  scattering  occurring  between  dif¬ 
fraction  peaks  Is  another  Indication  of  the  high  quality  of  the  overlayers. 

The  experimentally  derived  lattice  constants  for  each  of  the  overlayers  dis¬ 
cussed  In  this  paper  are  presented  In  Table  1.  In  general,  a  rare  gas  film 
could  be  used  for  several  hours  before  being  discarded.  When  each  experi¬ 
ment  was  completed,  the  coverage  was  checked  using  temperature  programmed 
thermal  desorption.  Since  the  monolayer  desorption  peak  always  occurs  at  a 
higher  temperature  than  that  of  the  multilayers  (due  to  the  magnitude  of  the 
rare  gas-Ag  holding  potential),  the  ratio  of  the  total  Integrated  area  to 
the  monolayer  peak  area  was  used  to  quantitatively  measure  the  coverage.  The 
TOF  data  collected  were  used  only  when  the  coverage  determined  fell  within  101 
of  that  expected. 

Oata  collection  was  performed  with  a  computer  controlled  multichannel 
scalar.  An  Incident  angle  of  45°  was  used  for  most  of  the  experiments 
reported  in  this  paper..  The  entire  surface  Brfllouln  zone  could  be  probed 
by  varying  the  detector  angle  while  holding  the  angle  of  Incidence  fixed. 
However,  phonon  energies  under  I  meV  could  not  be  well  resolved  due  to  inter¬ 
ference  from  the  elastically  scattered  peak.  Most  of  the  TOF  spectra  had 
to  be  signal  averaged  for  I  hour  In  order  to  achieve  a  reasonable  slanal-to- 
nolse  ratio.  This  was  necessitated  by  both  the  low  probability  of  inelastic 
scattering,  and  by  the  high  mass  4  background  in  the  ion-oumoed  detector. 


A  non-linear  least  squares  routine  was  used  to  accurately  fit  the  TOP 
spectra.  Phonon  energies  were  determined  by  comparing  the  peak  positions 
of  the  inelastic  transitions  to  the  time-of-fl  Ight  of  the  specular  beam. 

Although  using  peak  positions  Is  not  rigorously  correct,  the  error  intro¬ 
duced  by  this  procedure  Is  much  smaller  than  the  experimental  uncertainty, 
due  to  the  extremely  narrow  energy  spread  of  the  Incident  He  beam,  ^200ueV. 
Figure  2  presents  some  typical  TOF  spectra  and  fits  obtained  In  this  study. 

The  well  deffned  peaks  Indicate  that  single  phonon  exchanges  are  the  domin¬ 
ant  Inelastic  scattering  channel.  Each  TOF  spectrum  also  Includes  an  elastic 
feature,  which  we  attribute  to  Incoherent  scattering  from  crystal  defects. 

The  computer  simulations  In  Fig.  2  are  forward  convolutions  which  have  been 
scaled  to  show  excitation  llneshapes,  but  not  probabilities,  and  which  1) 
use  as  Input  the  experimental  dispersion  curves  11)  Include  the  Instrument 
transfer  function  and  11.1)  assume  single  phonon  exchanges  having  delta- 
function  energy  widths. 

Finally,  surface  phonon  dispersion  curves  were  constructed  by  plotting 
the  observed  phonon  energies  versus  parallel  momentum.  This  can  be  accomplished 
In  a  very  straightforward  way  using  the  conservation  equations  for  total  energy 
and  crystal  momentum.  Figure  3  shows  the  results  as  observed  along  the  fW 
direction  for  Ar(lll)  films  of  increasing  thickness.  The  dispersion  curves 
shown  In  this  figure  clearly  reveal  how  the  surface  dynamics  of  a  thin  rare 
gas  film  evolve,  on  a  layer-by-layer  basis,  from  monolayer  to  bulk.  The  mono- 
layer  mode  Is  dispersionless,  (characteristic  of  an  Einstein  oscillator  mode) 
while  the  bulk  (25  layers)  mode  will  be  shown  to  be  the  surface  Rayleigh  wave. 
Similar  dispersion  curves  have  been  measured  for  Kr  and  Xe.  The  comparison  of 
the  curves  to  the  results  of  lattice  dynamics  calculations  using  various  force 
laws  and  structures  as  Input  forms  the  central  theme  of  this  paper. 
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1 1 1.  Tha  Lattice  Dynamics  of  Saga  Gas  Overlavscs 

Th«  theory  o£  lattice  dynamic's  can  calculata  t ha  dispersion 
curves  of  phonons  by  solving  fot  fcha  eigenvalues  of  a  csrtaln 
matrix ,  known  as  tha  dynamical  matrix#  at  salactad  points 
along  tha  lina  of  intarast  in  tha  two  dimansional  surfaca 
Brillouin  zona.  Tha  eigenvalues  at  tha  point  ara  simply  tha 
sqoaras  of  tha  phonon  frequencies.  In  ordac  to  use  tha  tech¬ 
niques  of  lattica  dynamics,  thara  is  a  cartain  minimum  sat  of 

information  needed,  before  calculations  can  ba  dona.  This 

* 

information  falls  into  two  claasaa:  tha  first  is  tha  gaomatrical 
structura  of  tha  systam  (i.a.  tha  positions  of  tha  atoms),  and 
tha  sacond  daals  with  tha  intaraetions  batwaan  tha  atoms.  Given 
this  information,  than  lattica  dynamical  calculations  of  varying 
dagraas  or  sopmstication  can  oa  cons. 

For  tha  rare  gases  on. tha  Ag(lll)  surfaca,  tha  gasas  adsorb 
to  form  ordarad  multilayars,  as  mantionad  aarliar.  Thasa  layars 
hava  haxagonal  structura  parallel  to,  and  ara  incommensurate 
with,  tha  kg(lLl)  surface.  Thus,  for  tha  structural  inputs  to 
tha  lattica  dynamics,  wa  will  use  a  sat  of  layers  with  a 
haxagonal  lattica  structure.  For  tha  tri layer  case  wa  will 
consider  tha  affects  of  fee  stacking  of  tha  layars  (abc)  versus 
hep  stacking  (aba)  on  tha  dispersion  curves. 

Tha  kg  surfaca  appears  to  provide  primarily  a  support  on 


which  tha  rara  gas  crystal  can  ba  grown.  Tha  rara  gas-substrata 

interaction  is  a  potential  with  strong  dependence  on  the  height 

»  • 

of  tha  atom  above  tha  surface,  and  very  weak,  so  far  as  one  can 
sea  from  earlier  experiments,  dependence  on  the  atom's  lateral 


position.  The  a tom-subs.tr a ts  potential,  V(z)  ,  can  be  modeled  as 

a  van  der  Weals  1/z4  tail  at  large  z  combined  with  an  exponential 

repulsive  term  at  small  2.  This  is  the  form  that  Vidal i »  Cole 
12 

and  Klein  have  used  to  derive  a  proposed  universal  physisorp- 
tion  potential,  following  their  work,  we  take 

▼(*)  -  0  g(z#)  (1) 

where  z*  •  (z-za}/x.  Sere  g(o)  *  -  1,  so  that  0  is  the  depth  of 
the  potential  at  the  minimum,  zm  is  the  loeation  of  the  minimum 
of  the  potential  (which  is  the  position  where  the  atoms  in  the 
monolayer  sit) ,  and  l  is  a  parameter  with  units  of  length  defined 
by  l  •  (Cj/D) 1//3.  Here  C3  is  the  coefficient  of  the  1/z^  term  in 
V(x)  i.e.  V(x)  *  -  C3/z3  as  z  ♦  •  (We  have  tak'n  the  positive  z 
axis  to  the  normal  to  the  surface  and  pointing  into  the 
vacuum)  •  Tne  tunctionai.  tom  or  gtz#j  is  given  as 

*<*>  *  (&)  •'ax/'  *  — l-j  t2> 

The  parameters  u  and  a  are  dimensionless  and  are  related  by  the 
condition  that  g(8)  »  -1. 

In  order  to  model  the  atom-substrate  interaction,  . it  is 
necessary  to  determine  the  values  of  0,  C3,  1,  a,  and  u. 

Although  there  are  five  parameters,  there  are  only  three 
independent  ones  0,  C3,  and  a  (or  u) .  0  has  been  determined 

experimentally  by  Onguris,  at.  al.13  for  the  systems  of  interest 
here.  We  have  chosen  to  use  the  value  of  C3  calculated  by 
Zaremba  and  (Cohn.14  The  third  piece  of  information  that  we  have 
is  the  frequency  of  the  monolayer  vibrations  associated  with 
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adsorbate  motion  normal  to  the  surface.  This  is  provided  by  the 
data  reported  in  the  present  paper.  This  tells  us  the  second 
derivative  of  V(z)  at  its  minimum.  This,  together  with  0  and  C3 
is  sufficient  to  determine  a  and  hence  u.  The  V(z)  has  been 
determined,  and  we  can  then  incorporate  the  influence  of  V(z)  on 
the  lattice  dynamics.  For  the  three  adsorbates  of  interest. 

Table  1  provides  the  parameters  we  have  used. 

With  V ( z)  now  determined,  we  have  a  model  of  the  rare  gas 
atoms'  interaction  with  the  substrate,  but  before  any  calcula¬ 
tions  can  be  done,  we  need  to  know  how  the  rare  gas  atoms 
interact  with  each  other.  If  we  simply  wanted  to  model  the 
lattice  dynamics  of  the  system,  we  could  take  a  set  of  force 
constants  between  the  rare  gas  atoms,  and  determine  their  values 
by  adjusting  them  until  the  measured  dispersion  curves  have  been 
matched.  Rather  than  do  that,  we  have  chosen  to  make  use  of  some 
of  the  excellent  gas  phase  pair  potentials  that  have  been  deter¬ 
mined  by  Barker  and  co-workers.  The  potential  used  for  argon  was 
that  described  by  Barker,  Fisher  and  Watts .1S  The  potential  used 
for  krypton  was  the  Barker  K2  potential,16  and  the  Barker  X4 
potential16  was  used  for  xenon.  These  potentials  are  multipara¬ 
meter  in  nature,  and  reproduce  a  wide  range  of  physical  processes 
quite  well.  Since  they  have  analytical  forms,  it  is  quite  easy 
to  calculate  the  derivatives  required  for  the  lattice  dynamics. 
Although  these  potentials  reproduce  liquid  and  solid  properties 
as  well  as  gas  phase  properties  well,  there  should  be  some 
changes  in  the  pair  potentials  because  of  the  proximity  of  the 
adsorbates  to  the  Ag  surface.  However,  at  the  current  time  there 
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is  no  clear  experimental  information  available  on  what  these 

changes  might  be  and  available  theories  are  quite  approximate  in 

nature.  Hence  we  will  use  the  unmodified  pair  potentials  and  we 

may  then  regard  the  comparison  between  our  calculations  and  the 

data  as  a  test  of  whether  these  forms  may  be  applied  to  physi- 

sorbed  rare  gas  atoms.  This  analysis  then  provides  a  base  upon 

which  a  systematic  study  of  deviations  may  be  based.  Deviations 

from  the  gas  phase  potentials  may  be  produced  by  the  inclusion  of 

three  body  forces  such  as  Axilrod-Teller  forces,  and  substrate 

mediated  forces  such  as  those  of  McLachlan.^7  In  our  view,  among 

the  most  important  is  incorporation  of  the  influence  motion  of 

substrate  atoms  on  the  overlayer  lattice  dynamics  in  an  explicit 

fashion.  We  presently  have  this  problem  under  study,  and  we 

10 

describe  our  results  in  a  subsequent  publication. 

With  the  structure,  adsorbate-substrate  and  lateral  interac¬ 
tions  between  adsorbates  determined  as  just  outlined,  it  is  now 
possible  to  do  lattice  dynamical  calculations.  We  have  chosen  to 
explore  three  basic  models,  arranged  in  order  of  increasing 
sophistication.  The  first  model,  model  A,  assumes  that  only 
nearest  neighbor  central  force  interactions  exist.  Thus,  each 
atom  is  connected  by  "springs"  only  to  its  nearest  neighbors. 

The  nearest  neighbor  spacing  is  that  of  a  pair  of  atoms  in 
equilibrium.  The  only  coupling  between  the  substrate  and  the 
rare  gas  lattice  occurs  from  the  first  layer  (the  one  nearest  the 
surface)  to  the  substrate  with  the  spring  constant  ks  determined 
from  the  second  derivative  of  V(z).  The  rare  gas-rare  gas  force 
constant  k0  is  determined  from  the  second  derivative  of  the  pair 
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potentials ,  evaluated  at  the  minimum  of  the  potential .  This  is 
the  simplest  model  used. 

The  second  model#  model  3#  again  uses  only  neaseet  neighbor 

•  • 

central  forces  between  the  rare  gam  atoms.  Here  again  the 
nearest  neighbor  spacing  is  the  same  as  in  model  A.  The 
difference  is  that  now  every  layer  in  the  rare  gas  crystal  is 

et 

connected  to  the  substrate  by  a  force  constant  k_(jt  )#  where  i 

s  s  z 

labels  the  layers.  In  model  A#  k#U^  •  o  if  lz*  1.  Model  B 
takes  into  account  the  fact  that  the  atoms  in  the  higher  layers 
feel  V(s)j  we  compute  k gUa)  from  7(a). 

In  model  C  we  allow  all  the  rare  gas  layers  to  talas  both 
parallel  and  normal  to  the  silver  substrate#  to  assume  a 
configuration  in  which  the  static  potential  energy  is  minimised. 
At  the  same  time#  we  also  estsnd  the  effects  of  the  Barker 
potentials  out  beyond  the  nearest  neighbors.  This  is  accom¬ 
plished  by  generating  the  dynamical  natsis  by  a  computer  routine 
rather  than  setting  it  up  analytically.  Thus  any  number  of 
neighbors  can  be  used.  The  model  allows  us  to  see  the  effects  of 
the  long  range  tall  of  the  rase  gas  pair  potential#  and  also  of 

t 

the  selasation  effects  just  described.  The  number  of  neighbors 
used  is  shown  in  a  subscript  to  the  model#  in  the  subeequent 
discussion.  That  is#  the  Cj  model  which  includes  interactions 
through  fifth  neighbors#  and  so  on 

So  Figures  4  to  I#  we  show  the  results  of  the  lattice 
dynamics  calculations  for  the  three  rare  gases  on  Ag(111).  Parts  a,  b,  c 
and  d  show  the  dispersion  curves  for  the  monolayer,  bilayer,  trilayer  using 
Model  8,  and  2S  layers  using  Model  A,  respectively.  For  25  layers.  Model  A 
..in  nnt  Hif-fer  aooreclablv  from  Model  8.  The 
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curves  <:«  qualitatively  similar  foe  all  three  models.  We  shall 
elaborate  on  the  quantitative  differences  below.  For  the  mono* 
layer#  the  notion  normal  to  the  surface  is  decoupled  from  the 
notion  parallel  to  the  surface  in*  any  aodel  within  whieh  lateral 
interactions  between  adsorbates  have  central  fosee  character. 

This  property  is  shared  by  all  three  models.  Also#  the  two  nodes 
that  involve  motion  parallel  to  the  surface  decouple  into  a  pure 
transverse  node  and  a  pure  longitudinal  node  along  the  entire 
line  f  to  8.  Note  that  the  parallel  modes  are  acoustic  in 
nature#  while  the  motion  normal  to  the  surface  acts  like  an 
Sins tain  oscillator#  with  zero  dispersion. 

for  the  bilayer#  there  are  sir  modes  for  each  value  of  ^ • 

Two  of  these  are  transverse  modes  and  are  represented  by  the 
dashed  lines  is  the  figures.  The  other  four  nodes  are  sagittaiJ.y 
polarized  in  the  y-z  plane#  and  their  polarization  changes  as  one 
moves  along  the  f  -  3  line.  Let  os  label  the  modes  by  their 

polarization  at  the  f  point.  Then  these  are  two  longitudinal 

.  • 

modes  L^  and  Lj#  and  two  perpendicular  modes  and  i2#  la 

addition  to  the  two  transverse  nodes  T^  and  Tj-  At  f#  T^  and  L^ 

are  zero  frequency  acoustic  nodes  where  an  atom  in  each  layer 

% 

displaces  by  the  same  amount  is  the  appropriate  direction. 

At  f#  the  nodes  Tj  and  have  optical  character#  and  thus  the 
atoms  in  the  different  layers  neve  by  the  same  amount#  but  in 
opposite  directions.  The  node  has  atoms  in  the  different 
layers  all  noving  in  the  ♦  z  direction#  but  by  different  amounts, 
while  l2  has  the  atoms  noving  in  opposite  directions  130°  out  of 
phase#  again  with  different  magnitudes  of  displacement. 
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.  AC  8,  the  situation  is  lass  simple,  as.  the  atoms  in  the 
different  layers  no  longer  move  exactly  in  phase  or  out  of 
phase.  Here  T^  and  T2  &eoain  transverse  modes  with  the  atoms 
having  the  same  magnitude  displacements,  but  now  the  atoms  in  the 
second  layer  are  120°  in  phase  ahead  of  the  first  layer  atoms  in 
T2,  while  they  are  60°  behind  for  T2.  The  mode  is  now  no 
longer  a  simple  longitudinal  mode.  It  is  dominated  by  the  motion 
of  the  second  layer  atoms  in  the  z  direction.  In  contrast  to 
this*  L2  is  still  longitudinal,  and  is  dominated  by  the 

m 

longitudinal  motion  of  the  first  layer  atoms  along  the  y 
direction.  The  model  is  now  dominated  by  the  motion  of  the 
first  layer  atoms  normal  to  the  surface,  while  i2  is  dominated  by 
the  longitudinal  motion  of  the  second  layer  atoms.  For  all  of 
the  sagittal  modes,  the  other  layer  has  some  motion  in  the  y-z* 
plane,  but.  over  80%  of  the  contribution  to  the  normalization 
sum  l  |e  (ia)|*  comes  from  the  dominant  motion.  Here  e tt  2)  is 

V 

the  displacement  of  the  atom  la  layer  xs,  and  the  sum  just  quoted 
equals  unity.  Both' and  i2  have  the  second  layer  leading  the 
first  -by  120°,  while  b2  and  have  it  trailing  the  first  larger 
by  60°. 

The  situation  is  different  for  the  trilayer.  Here  the 
motions  are  exactly  in  phase  or  out  of  phases  at  f  and  also 
at  8  as  well.  At  F,  the  transverse  and  longitudinal  modes  form 
three  degenerate  pairs.  The  first  pair,  T^  and  L^,  are  uniform 
displacements  of  the  atoms  in  all  three  lyers  and  represent  rigid 
body  translations  of  the  rare  gas  crystal.  These  thus  have  zero 
frequency,  by  virtue  of  the  independence  of  V(z)  on  the  adsorbate 
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co-ordinates  in  the  plane  parallel  to  the  .surface.  The  second 
pair#  C.J  ar^  t2  a:a  analogous  to  the  symmetric  stretching  node  of 
a  linear  triatoaic  homo nuclear  molecule.  The  atoms  in  the  first 
and  third  layers  move  with  equal  magnitudes  and  in  opposite 
directions ,  while  the  atoae  in  the  second  layer  are  notionless* 
The  last  pair#  T3  and  C3,  are  analogs  to  the  asyaaetric  stretch, 
as  the  atoae  in  the  first  and  third  layers  nove  by  the  saae 
aaeunt  in  the  same  direction,  while  the  second  layer  atoae  move 
in  the  opposite  direction  with  twice  the  displaceaent*  The  nodes 
i^,  i 2*  4n^  I3  ell  have  the  atoae  aoving  along  the  z-direetion, 
with  having  all  the  atoae  noviog  in  the  saae  direction  by 
different  mounts*  At  fl  T^,  Tj,  and  T3  have  exactly  the  saxse 
eigenvectors  that  they  had  at  F.  Sowever,  is  now  doainated  by 
the  aoraal  notion  of  the  third  layer  atoae,  while  the  najor 
contribution  to  t.2  coaee  froa  the  noraal  notion  of  the  first 
layer  atoae*  Finally,  C13  has  the  eigenvector  that  L2  had 
at  f,  i*e*  it  is  now  the  analog  to  the  syaaetric  stretch, 
while  I3  is  now  doainated  by  the  longitudinal  notion  of  the 
second  layer  atoae* 

As  we  go  froa  3  layers  to  23  layers,  we  approach  the  case  of 
seai-infinite  argon.  How  we  see  a  single  Rayleigh  wave,  above 
which  the  bulk  bands  are  prtsant.  Note  the 

gap  in  the  bulk  bands  in  the  vicinity  of  3.  A  surface  mode  of 
primarily  longitudinal  polarization  is  present  in  the  gap.  Hints 
of  the  gap  can  be  seen  froa  the  dispersion  curves  for  the 
trilayer  case* 


Since  the  polarization  changes  as  on#  nova*  across  tho  zona, 
different  branches  can  ba  zasponsibla  for  axciting  tha 
experimentally  obsarvad  paaica  in  diffarant  parts  of  tha  zona, 
sinca  it  is  tha  perpendicular  atomic  motions  that  dominate  tha 
loss  cross  saction. 

Zn  Figs.  4  through  6,  .tha  chainad  linas  raprasant  a  bast  fit 
to  tha  azpariaantal  data.  . It  should  ba  remarked.  that  nodal  3 
doas  not  provide  ona  bast  fit  to  tha  data,  but  this  coapaxison 
will  prowl da  an  initial  oriaatation.  At  r,  tha  ix  is  rasponsibla 
for  tha  peak,  but  as  ona  moves  along  tha  lina  f  -  8,  a 
hybridisation  occurs  at  about  1/S  of  tha  way  to  8,  and  tha 
branch  becomes  tha  soda  axeitad  by  tha  atoms  from  thara  to  8. 

For  tha  As  and  Ke  bi layers,  thara  is  experimental  evidence  for  a 
second  paak  at  a  higher  frequency.  Ona  should  associate  thik 
peak  with  either  tha  lj^  or  tha  L?  branch;  we  shall  discuss  tha 
intaspratation  of  this  feature  later. 

Hawing  seen  something  of  tha  nature  of  tha  different  phonon 
branches,  and  how  their  character  varies  over  the  lina  f  -  8,  it 
is  useful  to  examine  tha  sensitivity  of  tha  theoretical  results 
to  tha  input  parameters  and  to  tha  choice  of  nodal.  We  present 
results  of  our  comparison  of  various  nodels  in  tha  fora  of  a 
sequence  of  tables  which  list  frequencies  at  F  and  8,  because 
soma  differences  are  too  small  to  show  up  in  dispersion  curve 
plots.  Table  3  shews  the  force  constants  (in  THz2)  used  for 
nodal  A,  including  two  choices  fac  k0,  one  from  tha. accurate  pair 
potentials  and  ona  from  tha  Cennard-Jones  6-12  potential,  for  tha 
thv*«  ?sr*  oases.  The  value  for  k«  was  defined  from  tha  experi- 
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mentally  measured  monolayer  frequency  via  kg  *  ^expt'  an<*  ot^er 
force  constants  defined  in  this  paper  are  in  equivalent  units. 
That  is,  in  the  case  where  the  force  constant  is  proportional  to 
the  second  derivative  of  the  two  particle  potential,  our  tabu* 
lated  force  constants  are  twice  this  second  derivative  divided  by 
the  adsorbate  mass.  It  is  immediately  obvious  that  the  Barker 
pair  potential  value  for  kQ  is  about  50%  larger  than  the  Lennard- 
Jones  value  for  kQ.  Hence  those  frequencies  which  depend  only  on 
kQ  should  be  about  20%  larger  for  the  Barker  pair  potential  then 
they  are  for  the  Lennard-Jones  potential.  This  can  be  seen  in 
Tables  3  through  8  by  comparing  the  columns  labelled  A  (B)  and 
A (LJ) •  Column  A(B)  gives  the  results  using  model  A  with  k0  from 
the  accurate  pair  potential,  while  A(LJ)  is  the  results  using 
model  A  with  the  Lennard-Jones  value  for  k0.  In  all  cases,  the 
frequencies  are  lower  in  the  column  A  (LJ)  than  those  in  column  A 
for  the  same  branch,  and  the  Lennard-Jones  results  are  in  clear 
disagreement  with  experiment. 

Another  interesting  issue  is  the  lateral  position  of  the  top 
layer  in  the  trilayer.  Does  it  sit  over  the  hep  or  fee  site? 

Bulk  rare  gas  crystals  form  fee  lattices,18  but  it  might  be 
possible  for  the  layer  to  sit  in  the  hep  site  for  the  tri layer 
system,  particularly  when  one  considers  that  energy  differences 
between  the  bulk  hep  and  bulk  fee  structure  are  very  small.  The 
results  for  the  hep  site  calculations  can  be  seen  in  Table  6  to  8 
under  the  column  A  (hep),  and  also  that  labeled  C5(hcp).  Notice 
that  at  f  and  R  there  is  no  difference  between  the  frequencies 
that  lattice  dynamics  predicts  for  the  fee  and  the  hep  models, 
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using  caodel  A,  and  only  small  changes  in  aodel  C.  As  one  novas 
away  from  these  high  sytnmatry  points  toward  tha  cantas  of'  tha  lina 
P  -  3#  than  soma  small  diffasaneas  appaar*,  but  tha  diffaraneas  ara 
so  small  that  tha  aspasimant  cannot  resolve  tham  at  peasant.  Thus 
m  can  oa km  no  dacision  on  this  quastion  at  the  present  tint. 

•  Whan  tha  atom-surface  potential  is  allowed  to  influence  tha 
atoms  in  tha  highes  layers ,  and  not  only  that  closest  to  tha 
surface/  it  creates  a  driving  force  to  causa  relaxation  between 
the  rare-gas  overlayers.  This  driving  fores  moves  the  higher 
planes  inward/  away  from  their  perfect  crystal  positions.  If  k0 
is  held  fixed/  then  the  resulting  changes  to  the  dispersion 
curves  are  small  and  occur  mainly  in  the  transverse  modes/  which 
have  not  yet  been  seen  experimentally.  More  important  are  the 
changes*  in  the  dynamical  matrix  for  the  higher  layers  caused  by 
7(2).  The  one  parameter  which  is  undetermined  in  the  Cole  and 
Vidal  1  form  of  7(2)  is  xB/  the  location  of  the  potential 

minimum.  This  presents  no  problem  for  the  monolayer/  .as  it  must 

•  • 

sit  at  Xg.  Hence  we  can  choose  the  minimum  to  be  the  origin 
(2«0)  and  the  value  of  2,  is  irrelevant.  Siut  foe  multilayer 
structures/  we  must  allow  the  layers  to  relax/  and  the  innermost 
layer  no  longer  sits  at  a  minimum  of  7(z).  We  have  performed  a 
simple  calculation  to  get  a  rough  idea  of  the  amount  of  inter¬ 
layer  relaxation  that  occurs.  We  have  calculated  the  energy  per 
unit  cell  by  summing  the  pair  potentials  over  shells  of  atoms 
until  the  sum  converges/  and  then  added  the  contribqtion  from 


7(s) •  We  varied  the  two  dimensional  lattice  constant  a0  and  the 
location  of  the  layers  to  locate  the  minimum  in  the  total 


energy.  The  results  arc  summarised  in  Tabic  10.  As  sight  be 
expected,  the  first  layer  is  pushed  closer  to  the  substrate  by 
the  layers  above  it.  Knowing  the  position  of  the  layers,  it  is 
then  possible  to  calculate  the  force  constant  acting  on  that 
layer  from  7(a)  •  The  ksU2)  «•  shown  in  Table  11.  Tt*eir  effec 
can  be  seen  in  Tables  4  through  9  under  the  column  labelled  a. 
Sote  that  in  model  B  the  interplanar  spacings  are  not  changed, 
only  the  forces  connecting  the  planes  to  the  substrate  are 
altered.  It  can  be  seen  that  only  those  modes  with  vibration 
normal  to  the  surface  are  affected,  when  these  results  are 
compared  with  column  The  degree  to  which  they  are  affected 
depends  directly  on  the  amount  of  normal  vibration  they  have, 
hence  the  and  1^  and  fl  modes  are  the  most  strongly  affected. 

Since  the  rare  gas  pair  potentials  are  fairly  long  rangeo, 
one  must  inquire  about  the  role  of  distant  neighbor  interactions 
We  thus  turn  to  a  discussion  of  the  calculations  based  on  model 

C.  It  was  found  that  the  convergence  was  very  rapid,  as  we 

•  • 

included  progressively  more  neighbors.  The  most  important 
change,  however,  was  the  fact  that  this  method  incorporates 
changes  in  the  nearest  neighbor  distances.  * In  the  nearest 
neighbor  models  A  and  B,  the  atoms  within  a  given  plane  have  as 
their  equilibrium  position  the  minimum  of  the  pair  potential. 
Hence  Jc0  is  calculated  from  the  second  derivative  of  the  pais 
potential  evaluated  at  its  minimum.  When  we  use  model  C,  the 
atoms  are  no  longer  at  the  minimum.  They  are  at  the#  locations 
shown  in  Table  10.  The  lattice  constant  is  smaller  than  the 
distance  at  which  the  pair  potential  has  its  minimum,  since  the 


lattice  can  lower  its  total  energy  this  way.  Table  12  shows  the  difference 
that  this  makes  for  the  A r  bi layer  by  going  from  model  3,  to  -odel  C-j  (near¬ 
est  neighbors,  but  with  a  new  kQ).  We  see  that  the  frequencies  are  raised 
(with  the  increased  kQ)  at  the  reduced  separation  of  the  atoms.  As  we  add 
more  neighbors  to  the  calculation  of  model  C  we  see  the  frequencies  are  now 
reduced  from  the  values  given  by  model  C^.  Note  that  convergence  Is  rapid, 
and  that  by  the  time  five  neighbors  have  been  used,  the  frequencies  have 
converged  to  four  figures.  Finally,  we  present  in  Table  13  the  normal  mode 
frequencies  at  1?  for  the  Ar,  Kr,  and  Xe  Rayleigh  waves  calculated  using  models 
A  and  C^. 

The  last  general  consideration  we  discuss  Is  the  effect  of  relaxation 
normal  to  the  substrate  of  the  multi -layer  structure,  within  the  framework  of 
model  C.  In  model  C,  the  force  constants  will  vary  as  the  crystal  relaxes,  and 
thlshas  the  possibility  of  Influencing  the  frequencies  considerably.  This  can 
be  seen  by  comparing  the  relaxed  results  (column  Cg  (r))  and  the  unrelaxed 
results  (colimm  Cg)  In  Tables  4  to  9.  Changes  in  frequency  as  large  as  ten 
percent  can  occur  as  a  consequence  of  interlayer  relaxation. 

IV.  Discussion  and  Concluding  Remarks 

A  result  of  our  comparison  between  the  various  models  is  that  the 

■  t 

overall  structure  of  the  dispersion  curves  is  rather  insensitive  to 
the  manner  In  which  the  lateral  interactions  between  adsorbates  is  modeled; 
clearly  use  of  the  lennard-Oones  potential  gives  results  substantially  at 
variance  with  the  other  models  explored,  however.  By  and  large,  a  nearest 
neighbor  central  force  model,  possibly  with  force  constant  calculated 
from  the  potential  of  Barker  and  co-workers,  provides  a  description  of 
the  dispersion  curves  fully  adequate  to  account  for  the  data 


presently  available.  In  fact  modal  C5#  which  gives  the  best 
agreement  with  experiment,  yields  frequencies  at  M  and  f  which 
ere  at  most  6%  higher  than  frequencies  o£  model  A  for  those  modes 
which  are  detected  in  the  experiments.  Relaxation  effects  also 
have  a  modest  influence  on  the  phonon  spectra.  To  repeat  an 
earlier  comment#  we  believe  that  coupling  between  the  adsorbates# 
and  substrate  atom  motions#  is  potentially  more  .important  than 
modest  variations  in  the  nature  of  the  adatom-adatom  potential. 

In  our  view#  earlier  work  in  the  area  has  perhaps  over  emphasised 
this  last  issue#  though  unambiguous  determination  of  substrate- 
induced  modifications  to  the  lateral  interactions  between  physi- 
sorbed’ atoms  is  surely  an  issue  of  fundamental  interest,  me  now 
turn  to  a  more  detailed  discussion  of  the  comparison  between 
theory  and  experiment. 

The  most  sophisticated  model  used  above  is  model  C#  where  the 
phonon  dispersion  relations  are  calculated  by  constructing  a 
dynamical  matrix  about  a  static  configuration  which  minimizes  the 
total  energy  of  the  model  system.  When  the  predictions  of  this 
model  are  compared  with  the  data#  there  are  in  fact  some  signifi¬ 
cant  discrepancies  between  theory  and  experiment.  For  the  Ar 
bi layer#  the  calculated  width  of  the  phonon  band  is  distinctly 
too  large#  and  we  find  a  similar  problem  with  all  the  trilayers. 

The  lattice  constants  we  calculate  by  minimizing  the  total 
energy  of  model  C  are  in  fact  slightly  smaller  than  those 
measured  in  the  experiment.  This  may  be  seen  from  Table  10.  The 

•  i 

discrepancies  are  largest  for  the  Ar  bi layer#  and  all  the 
trilayers#  and  it  is  in  these  cases  we  have  the  greatest  diffi- 
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culty  fitting  the  phonon  dispersion  curves.  In  fact,  in  our 
model,  we  have  omitted  all  three  body  contributions  to  both  the 
static  energy,  and  to  the  dynamical  matrix  used  in  the  lattice 
dynamics.  We  have  verified  explicitly  that  the  phonon  fre¬ 
quencies  are  affected  only  very  slightly  by  the  inclusion  of 
three  body  interactions.  We  have,  for  example,  assessed  the 
contribution  of  the  Axilrod-Teller  triple  dipole  interaction,  to 
find  its  contribution  is  quite  negligible.  However,  from 
theoretical  studies  of  bulk  crystals,  it  is  established  that 
three  body  couplings  are  required  to  generate  the  proper 
lattice  constant.  Our  problem  here,  if  we  wish  to  fully  include 
them  in  the  analysis,  is  that  we  have  no  information  on  non¬ 
central  force  contributions  to  the  energy  that  are  unique  to  the 
adsorbed  state.  Me  envision  three  body  interactions,  in  which 
one  of  the  three  entities  is  a  substrate  atom. 

However,  if  it  is  indeed  true  that  the  phonon  frequencies  are 
influenced  only  very  modestly  by  the  three  body  couplings,  then 
we  may  calculate  phonon  dispersion  curves  by  simply  using  the 
measured  lattice  constant,  rather  than  by  using  that  provided  by 
our  energy  minimization  calculation.  As  seen  from  Table  10,  the 
experimental  lattice  constant  is  a  bit  larger  than  we  calculate 
with  our  procedure,  we  argue  that  a  small  expansion  of  the 
lattice  parallel  to  the  surface  is  driven  by  three  body  inter¬ 
actions,  and  other  lateral  interactions  ignored  in  our  model.  We 
have  performed  calculations  using  the  experimental  lattice 
constant,  but  with  interlayer  spacing  relaxed  in  the  manner 
provided  by  our  model;  this  picture  does  not  provide  an  optimum 
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account  of  tha  data.  Tha  bast  agraanant  is  found  by  using  tha 
axparimantal  lattica  constant,  but  than  choosing  an  intaclayar 
spacing  appropriata  to  a  pat fact  fee  lattica.  Tha  lack  of 
vartical  talaxation  has  its  likaly  origin  in  setaaning  of  tha 
holding  potantial  by  tha  innar  adsoebata  layacs,  tha  intarlayar 
spacing  may  ba  axpandad  a  bit  foe  tha  sama  raason  (a.g.  theaa 
body  intaractions)  as  tha  two  dimansional  unit  call  is  axpandad. 
It  is  possibla  that  a  snail  amount  of  intarlayar  ralaxation  is 
prasant;  tha  prasant  pracision  in  tha  data  paraits  no  claar 
conclusion  on  this  point.  Ibis  picture  provides  an  excellent  fit  to  the 
data.  In  Fig.  7,  we  show  a  comparison  between  theory  and  experiment,  for  the 


tbrtt  di layer  sy stains.  The  crosses  are  the  data  points,  the  solid  line  a 

Pest  fit  to  the  data,  and  the  triangles  are  provided  by  the  theory,  assuming 

the  and  modes  are  responsible  for  the  loss  feature,  as  described  earlier. 

Fig.  3  shows  a  similar  comparison  for  the  three  trflayers,  and  Fig.  9  for 

the  twenty-five  layer  slab.  We  also  show  In  Tables  U-16  (again  using  the 

2 

experimental  lattice  constants)  the  values  of  |ez|  for  the  various  modes  that 
2,  3,  and  2S-layer  films  of  A r,  K r,  and  Xa/Ag(lll)  have  along  the  F*ff  direction. 

That  meet  significant  discrepancy  between  theory  experi¬ 
ment  occurs  for  the  xenon  txilayer,  where  the  theory  eyeteneti- 
eally  falls  below  the  data,  in  the  vicinity  of  F.  We  have  no 
clear  suggestion  to  offer  about  the  origin  of  this  difference, 
which  amounts  to  roughly  0-2  me?  for  the  data  cloeeet  to  f .  We 
are  investigating  the  peeaibility  that  coupling  between  the 
snhetrate  and  adsorbate  motions  may  raise  the  frequency  of  the 
acsoroata  mooes  near  r . 

For  the  Ar  and  tr  bi layer,  a  week,  loae  faature  is  observed  at 
frequencies  considerably  higher  than  the  range  eoversd  by  Fig. 

7  •  The  structure  is  seen  near  the  3  point  for  both  cases.  In 
Fig.  10  ,  the  crosses  indicate  these  features,  and  on  the  plot 
we  superimposed  the  dispersion  relations  of  the  i^  and  Lj  modes, 
both  of  which  have  displacement  components  normal  to  the  surface. 


and  thus  can  be  expected  to  produce  loss  features  in  this 
spectral  range.  In  fact,  when  we  calculate  |e^|2,  the  square  of 
the  eigenvector  component  normal  to  the  surface  for  the  outermost 
adsorbate  layer,  the  two  are  comparable  in  magnitude  for  these 
two  modes.  Thus,  we  argue  the  loss  feature  is  produced  by  the 
combined  action  of  those  two  eigenmodes,  and  is  an  unresolved 
doublet.  If  we  construct  a  single  loss  feature  by  superimposing 


26 


two  rather  diffuse  Lorentzians ,  and  assume  the  integrated 
strength  of  each  sealu  Uka  tho  square  of  tho  eigenvector  jez|2, 
wo  ox poet  tho  pook  of  tho  resulting  structure  to  oeear  at  tho 
frequency  5  •  («l|o^l) J2  ♦  | 2)/2( !«^1} | 2  *  (e^2) j  2)  .  For 

tho  argon  and  kypton  bilayer,  tho  froqaoney  Z  at  8  ia  indicatod 
by  an  arrow  plaeod  on  tho  plot.  Tho  froqaoneioa  weighted  in  thia 
fashion  aro  in  very  good  accord  with  tho  position  of  tho  two  loss 
foaturos  at  ».  At  Q/^  -  0.5,  has  tho  largest  value  of  |e2|2. 

Za  our  view,  a  Tory  intriguing  aspect  of  tho  prosont  work  is 
that  it  shows  vividly  how.  tho  Rayleigh  wave  of  tho  semi-infinite 
solid  evolves ,  as  a  function  of  layor  thickness,  out  of  tho 
nozaal  sodas  of  a  simple  monolayer.19’20  This  Is  Illustrated  very 
beautifully  ia  Figure  3  •  Tho  Sayloigh  wave  on  tho  semi- 
infinite  solid  is  always  an  adairture  of  motions  parallel  to,  and 
nozaal  to  tho  surface.  For  tho  monolayer,  with  lateral  interac¬ 
tions  of  central  force  character  only,  tho  normal  modes  aro 
polarized  either  parallel  to  tho  surface  (2  branches) ,  or  normal 
to  tho  surface,  when  coupling  to  substrate  motions  is  ignored. 
Thus,  since  Ho  atoms  scatter  predominantly  from  tho  normal 
component  of  displacement,  wo  see  tho  Tirtually  disperionless 
monolayer  mode  in  Fig.  3  .  For  multilayer  structures,  except  at 
high  symmetry  points,  tho  modes  polarized  in  the  sagittal  plane 
hate  mixed  character,  except  at  the  high  symetry  points.  The 
mode  that  evolves  into  the  Rayleigh  wave  acquires  dispersion,  and 
the  frequency  at  f  drops  monotonically  as  the  film  ..licknes. 
decreases,  to  zero  in  the  limit  of  infinite  thickness.  If  we 
consider  a  thick  film,  with  thickness  d,  then  the  dispersion 


relation  of  the  low  frequency  Rayleigh  aode  will  track  that 
appropriate  to  the  solid  in  the  regime  where  Q? d  »  1,  where 
is  the  wave  vector  of  the  wave,  parallel  to  the  surface.  In  this 
regiae,  the  wave  is  sufficiently  localised  to  the  adlayer/vacuua 
interface  that  it  does  not  "feel"  the  adlayer/substratq  inter* 
face.  Rhea  Qfd  becomes  the  order  of  unity,  then  the  wave 
penetrates  deeply  enough  to  feel  the  adlayer/substrate  interface, 
and  its  frequency  saturates  to  a  finite  value.  At  Qf  •  0,  we 

m 

have  a  pure  s  polarised  aode,  and  one  aay  argue  (for  a  fila  aany 
atoaic  layers  thick)  that  its  frequency  should  be  sc^/Zd,  where 
Cfc  is  the  appropriate  transverse  sound  velocity.  If  we  apply 
this  simple  formula  to  the  trilayer  and  the  bi layer,  then  as 
Qj  ♦  o,  the.  ratio  of  the  frequency  of  the  trilayer  aode  to  that 
of  the  bilayer  would  be  predicted  to  be  2/3,  wnicn  is  remarraoxy 
close  to  the  ratio  we  see  in  Fig.  3  . 

In  summary,  we  have  Investigated  the  surface  phonon  spectroscopy  of  rare 

gas  overlayers  adsorbed  on  Ag(111)  using  Inelastic  He  scattering  as  the  dynamical 

•  • 

probe.  This  was  done  on  a  layer-by-layer  basis  for  1,  2,  3,  and  25-layer  (111) 
oriented  films  of  Ar,  Kr,  and  Xe.  The  monolayers  exhibited  a  dispersion  less 
Einstein  oscillator  mode,  with  the  rare  gas  atoms  vibrating  In  an  uncorrelated 
manner  In  the  holding  potential  of  the  Ag  substrate.  As  progressively  more 
layers  were  added,  the  observed  modes  showed  Increasing  dispersion,  until  at  25 
layers  a  bulk-like  Rayleigh  mode  was  observed.  This  evolution  Is  due  to  the 
Increase  In  surface-to-substrate  distance  as  the  film  gets  progressively 
thicker,  making  the  Ag-rare  gas  interaction  less  Important  with  respect  to  the 
weaker  rare  gas-rare  gas  potential  In  determining  the  dynamical  behavior  of  the 


surface.  The  experimental  results  were  compared  with  lattice  cynamics  cal¬ 
culations  using  realistic  rare  gas  pair  potentials  derived  from  gas  phase 
experiments,  and  a  Ag-rare  gas  potential  derived  In  part  from  the  Observed  mono- 
layer  frequencies.  It  Is  important  to  note  that  only  the  laterally  averaged 
Ag-rare  gas  potential  Is  necessary,  since  all  of  these  rare  gas  systmns.  are 
translationally  Incommensurate  with  respect  to  the  Ag  substrate.  When  all 
layers  are  coupled  to  the  substrate,  and  the  rare  gas-rare  gas  interactions  are 
modelled  with  a  .realistic  gas  phase  pair  potential  summed  over  all  neighbors 
necessary  for  convergence,  very  good  agreement  Is  reached  between  theory  and 
experiment. 
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Experimental  lattice  constants,  aQ,  for  the  systems  stuci ad.  Results 
are  accurate  to  •  0.02  X. 

Table  2 

The  input  parameters  used  to  cone t rue t  V(z) ,  fes  the  three 
edeorbetee  of  iatereat  la  the  preaeat  paper. 

Table  3 

.  Input  parameters  for  model  A  oaed  to  deacribe  the  rare  gas 
over layers.  The  aeeoad  column  la  the  value  of  It,  provided  by  the 
Cole-Vidali  potential#  obtained  by  fitting  to  the  measured  value 
of  the  aoao layer  vibration  normal  to  the  aurfaee.  The  two  valuea 
of  k0  are  deduced  from  the  pair  potantiala  proposed  by  Barker  et. 

al.  (third  column) #  and  that  determined  from  the  Lennard-Jones 

•  .  » 

potential  (fourth  column) •  The  final  column  lists  valuea  of  k0 
obtained  empirically  by  fitting  the  measured  phonon  frequency 
at  5  tor  a  twenty  rive  layer  riim  to  coat  generatea  oy  moaei  a 
applied  to  such  a  slab. 

Table  4 

The  normal  node  frequencies  calculated  for  the  various  models 

discussed  in  the  teat#  at  f  and  H  in  the  two  dimensional  Brillouin 

rone  of  the  Ar  bilayer.  The  column  A(B)  are  frequencies 

calculated  for  model-  A#  with  the  foree  constant  between  rare  gas 

% 

atoms  determined  by  the  Barker  gas  phase  potential#  the  column 
A(LJ)  uses  the  Leonard  Jones  potential  to  generate  this  force 
constant  (see  Table  3) #  that  labeled  B  is  model  B#  and  two 
versions  of  model  C  are  eaplored.  Each  sums  over  pairs  out  to  a 
distance  equal  to  five  nearest  neighbor*,  the  first  (column 
labeled  Cj)  treats  an  unrelaxed  bilayer#  and  one  (column  labeled 
C5(r))  treats  a  fully  relaxed  system.  All  frequencies  are  in  mev. 
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Table  5 

The  normal  soda  frequencies  at  F  aad  A  la  the  two  dimensional 
Brillouia  zona,  for  the  Kr  bilayar.  Tha  columns  ara  labalad  by 
tha  convention  used  ia  Tab la  4 . 

Tabla  6 

. Tha  normal  mode  frequencies  at  F  aad  A  la  tha  two  diaaaaloaal 

Brillouin  zone ,  for  tha  Za  bilayar.  Tha  columns  ara  labalad  by 

*  • 

tha  convention  aaad  ia  Tabla  4. 

Tabla  7 

Tha  normal  aoda  fraquanciaa  at  F  aad  A  for  tha  Ar  trilayar, 
for  tha  thraa  aodala  explored  ia  this  work,  We  also  calculate  tha 
fraqoaaelaa  for  a  trilayar  with  HC9  stacking  arraagaaaat  ABAr 
rathar  than  tha  ICC  arraagaaaat  ABC.  This  has  baea  dona  for  nodal 
A  aad  nodal  Cj*  aad  taa  results  ror  aaea  ara  in  tna  columns 
labalad  A(HC?)  aad  C9(BC?X.  Baa  tha  caption  of  Tabla  4  fbr  tha 
r anaiaiag  conventions. 

Tabla  8 

Tha  normal  aoda  fraqoaaelaa  at  F  aad  A  for  tha  Kr  trilayar , 
for  tha  thraa  models  used.  Tha  columas  ara  labalad  by  tha 
convention  used  ia  Tabla  7. 

Tabla  9 

Tha  normal  mode  frequencies  at  F  and  A  for  tha  Xa  trilayar, 
for  tha  thraa  models  used.  Tha  columns  ara  labalad  by  tha 
convention  used  in  Tabla  7. 

Tabla  10 

Tha  affect  on  static  properties  of  tha  bilyars  and  tri layers 
of  extending  tha  influence  of  7(r)  beyond  tha  layer  closest  to 


the  substrate..  Here  in  the  first  column  z(l)  is  the  position  of  the  inner- 
most  layer,  measured  relative  to  z^,  d^  and  are  the  spacings  between 
the  first  and  second,  and  the  second  and  third  layers,  respectively,  and  aQ 
is  the  lattice  constant  calculated  for  the  structure.  The  last  column  gives 
the  experimentally  measured  lattice  constants. 

Table  11 

The  force  constant  ksUz)  for  the  various  layers,  generated  by  V(z) 
after  the  structures  are  fully  relaxed.  The  units  are  as  described  in  the 
text. 

Table  12 

For  the  argon  bilayer,  we  show  frequences  at  F  and  ft  calculated  with 
model  3,  and  with  various  versions  of  model  C.  The  subscript  on  C  in  the 
various  coitions  shows  the  number  of  neighbor  shells  Included  in  the  dynamical 
matrix. 

Table  13 

The  normal  mode  frequencies  at  ft  using  models  A(B)  and  Cg  for  the 
Rayleigh  wave  of  bulk  Ar,  Kr,  and  Xe. 

Table  14 

2 

Values  of  |ez|  ,  the  square  of  the  eigenvector  component  normal  to  the 
surface,  for  Ar  bilayer,  trilayer,  and  25-layer  films  (Rayleigh  wave  only). 

In  the  case  of  the  monolayer,  |e2j  is  equal  to  1  for  the  X  mode,  and  0  for 
the  other  modes. 

Table  15 

Same  as  Table  14,  but  for  Kr. 

Table  16 

Same  ^s  Table  14*  but  for  Xe. 
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Figure  1:  Diffraction  spectra  for  the  Ar  overlayers  studied. 

Experimental  results  ara  shown  for  that  3  -  45°  in  tha  (112) 
direction,  T#  •  21°  X,  incident  baas  energy  of  13  a»7,  and 
with  tha  background  substzacted.  Conditions  for  all  of  tha 
spectra  ara  tha  same,  so  intensities  ara  directly  comparable. 
Arrows  indicate  tha  angles  at  which  tha  phonon  spectra  of 
figure  2  ware  taken. 

floure  2:  TOT  spectra  of  tha  Ar  overlayers  studied,  Ts  •  21  K. 

and  incident  beam  energy  of  18  meV.  Tha  dots  are  tha 
experimental  data,  solid*  lines  are  a  least-squares  fit  of 
these  data,  and  dashed  lines  are  computer  simulations  (see  text). 
Arrows  Indicate  position  of  elastic  tlae  of  flight.  Spectra  were  taken  with 
1  usee  channels  and  one  hour  signal  averaglnq. 

Figure  3s  Dispersion  curves  of  the  Ar  overlayers  from  f  to  3. 
Open  figures  are  for  energy  loss  of  the  beam,  closed  figures 
are  for  energy  gain.  Representative  error  bars  are  shown  for 

i 

energy  and  momentum  loss  features  (gain  and  loss  are  with 
respect  to  the  incident  beam) ,  which  occured  in  the  f irst 
sone  (normal  processes).  Solid  lines  are  polynomial  fits  to 
the  data. 

Figure  4:  For  argon,  we  show  the  phonon  dispersion  curves 


calculation  from  model  B  for  (a)  the  monolayer,  (b)  the 
bilayer,  (c)  the  trilayer,  and  (d)  a  "bulk"  crystal.  The 


SI 


information  here  is  synthesized  from  a  calculation  based  on  a  twenty 
five  layer  slab.  For  l,  2,  and.  3  layers,  the  solid  lines  are  calcu¬ 
lated  dispersion  curves  for  modes  polarized  in  the  sagittal  plane,  the 
dashed  lines  are  transverse  modes,  and  the  chained  line  is  a  best  fit 
to  the  data.  For  25  layers,  triangles  are  calculated  dispersion  curves, 
and  the  chained  line  is  a  best  fit  to  the  data. 

Figure  5:  The  same  as  Figure  4,  but  for  krypton. 

Figure  6:  The  same  as  Figure  4,  but  for  xenon. 

Figure  7:  A  comparison  between  theory  and  experiment  fdr  the  three  bl layers. 
The  crosses  are  the  experlmntal  points.  The  slashes  through  selected 
crosses  are  indications  of  the  error  In  the  data,  lhe  solid  line  Is  a 
best  fit  to  the  data,  while  the  triangles  are  obtained  from  the  theory, 
assuming  the  modes  1*^  and  are  responsible  for  the  loss  peak,  as 
discussed  earller-ln  the  text. 

Figure  8:  The  same  as  Fig.  7,  but  fdr  the  sequence  of  trl layers. 

Figure  9:  The  same  as  Fig.  7,  but  for  a  twenty  five  layer  slab. 

Figure  10:  The  dispersion  relation  of  the  high  frequency  loss  features  near 
ff,  observed  for  both  the  argon  and  the  krypton  bl  layers.  Superimposed 
on  the  data  (.crosses)  are  the  dispersion  relations  of  the  and 

modes;  the  arrows  indicate  the  value  of  the  frequency  «  defined  in  the  text. 
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